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Abstract. Divide-and-conquer is a central paradigm for the design of algorithms, through which some 
fundamental computational problems, such as sorting arrays and computing convex hulls, are solved 
in optimal time within 0(n\ogri) in the worst case over instances of size n. A finer analysis of those 
problems yields complexities within 0(n(l + "H(ru,..., ru,))) C 0(n(l+ log k)) C O(nlogn) in the 
worst case over all instances of size n composed of k “easy” fragments of respective sizes n\,... ,rik 
summing to n, where the entropy function T-L{n \,. .., rik) = Yli=i yr 1°6 measures the “difficulty” 
of the instance. We consider whether such refined analysis can be applied to other algorithms based 
on divide-and-conquer, such as polynomial multiplication, input-order adaptive computation of convex 
hulls in 2D and 3D, and computation of Delaunay triangulations. 
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1 Introduction 


The divide-and-conquer paradigm is used to solve central computational problems such as Sorting ar- 
computing Convex Hulls 9;, Delaunay Triangulations and Voronoi Diagrams [6 of 

Polynomial Multiplica- 


rays 


10 


points in the plane and in higher dimensions, Matrix Multiplication 
tion 


19 


16 , Half-Plane Intersection [Hu , among others. For some of these problems, this paradigm yields 


an optimal running time within 0(nlogn) in the worst case over instances of n elements. 

An adaptive analysis of slight variants of some of these algorithms yields improved running times on large 
classes of instances. Those results can be refined up to complexities within 0{n{l+T-L{ni,... C 

0(n(l+ log fc)) C O(nlogn) in the worst case over all instances of size n composed of k “easy” fragments 
of respective sizes n,\, ..., n± summing to n , where the entropy function ..., n^) = Xu=i ~ l°g yy 

measures the “difficulty” of the instance. We describe here two examples for the problem of Sorting arrays 
(see Section [2] for more). 


— Munro and Spira [14| showed that the algorithm MergeSort can be adapted to sort a multiset S of n 
elements in time within 0(n(l+7i(mi, ..., rn a ))) C 0(n( 1+ logu)) C O(nlogn), where a is the number 
of distinct elements in S and mi,..., m a (such that J2i =1 m i = n ) are the p multiplicities of the distinct 
elements in S, respectively. 

— Taking advantage of the input order, Knuth [10 considered sequences formed by runs i.e., contiguous 
increasing subsequences, and described an algorithm sorting such sequences in time within 0(n(l + 
logp)) C O(nlgn), where p is the number of runs in the sequence. Barbay and Navarro (21 improved the 
analysis of this algorithm in time within 0(n( 1 + 'H(r 1 ,..., r p ))) C 0(n(l+ log p)) C O(nlogn), where 
ri,..., r p (such that Xu=i r i = n ) are the sizes of these runs. 
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Similar analysis techniques have been applied to some other problems, though only partially. Following an 


example for the problem of computing the Convex Hull in the plane (see Section 3.2 for more): Levcopoulos 
et al. 11 described an adaptive algorithm for computing the Convex Hull of a polygonal chain. The 


algorithm takes advantage of the minimum number k of simple subchains into which the polygonal chain can 
be partitioned. They showed that the time complexity of the algorithm is within 0(n(l+log k)) C O(nlogn). 

Hypothesis. Which similar refinements can be applied to which divide-and-conquer algorithms, if any, and 
to what depth? 

Our Results. In Section [2] we list and classify previous refined analyses between those that are Structure 
Based and those Input-Order Based. In Section 3.1 we describe a refined analysis of the principal step in the 
algorithm for Polynomial Multiplication using the Fast Fourier Transformation. In Sections |3.2| and 
|3.3[ we describe two distinct refined analyses for problems in computational geometry, which yield various 
optimal input-order adaptive results on the computation of Convex Hulls, Delaunay Triangulations, 
and Voronoi Diagrams in the plane. In Section 3.2 we refine the analysis of Levcopoulos et al. [Til’s 


algorithm for computing the Convex Hull of polygonal chains in time within 0(n( 1 + H(ni ,..., n K ))) C 
0(n(l+logre))) C O(nlogn), where m ,... ,n K are the lengths of the subchains of a partition of a polygonal 
chain of n points into the minimum number k of simple subchains. In Section |3.3[ we describe a refined 
analysis of the computation of Voronoi Diagrams and Delaunay Triangulations for sequences S 
formed by n points, which yields a time complexity within 0(n( 1 + TL[v\, ..., v M ))) C 0(n(l+ log p)) C 
O(nlogn), where ui,...,u M are the sizes of the minimum number p, of monotone histograms in which S 
can be cut, with respect to two fixed orthogonal directions and show, as a corollary, an upper bound for 
computing the CONVEX Hull of a sequence S formed by n points in time within 0{n(\ + %(v\, ..., v^))) C 
0(n(l+ log / T )) C 0(n log n ), where V \ ,..., are the sizes of the minimum number of monotone histograms 
in which S can be cut, with respect to two fixed orthogonal directions. In Section [4j we describe some more 
difficult applications of such refined analyses. In Section [5j we discuss the possibility of designing algorithms 
which analyses combine Structure Based and Input-Order Based results synergically as opposed to running 
them in parallel. 


2 Classification of Results 


We review here some results on the refined analysis of algorithms for Sorting arrays and for computing 
planar Convex Hulls, Delaunay Triangulations and Voronoi Diagrams. We classify the various 
refined analysis between those focusing on the structure of the instance (Section [2. 1[ ) versus those focusing 
on the order in which the input is given (Section 2.21. There does not seem to be any example where both 
strategies are mixed: we discuss the potential for those in Section [5] 


2.1 Structure Based Results 

By “Structure Based Results” we mean algorithms taking advantage of the structure of the instance, for 
example, taking advantage of the frequencies of the values in a multiset or of the relative positions of the 
points in a set. Such results are known for Sorting multisets and for computing Convex Hulls. 

MergeSort is a divide-and-conquer Sorting algorithm in the comparison model. This algorithm relies 
on a linear time merge process, that combines two ordered sequences into a single ordered sequence. 


Concerning the problem of Sorting, Munro and Spira 14 considered the task of Sorting a multiset 
S = {xi ,..., x n } of n real numbers with er distinct values, of multiplicities mi,... ,m CT , respectively, so 
that m i = n - They showed that adding counters to various classical algorithms (among which the 

divide-and-conquer based algorithm MergeSort) yields a time complexity within ..., to ct ))) C 

0(n(l+\oga)) C O(nlogn) for Sorting a multiset, where 7t{mi,... ,m CT ) = Yli = 1 — log — measures the 
entropy of the distribution of the multiplicities (mi,... ,m a ). This result takes advantage of the frequencies 
of the values i.e., the structure of the instance. 








Given a set P of n points, the Convex Hull of P is the smallest convex set containing P (see Figure [T]?? 
Considering the problem of computing the Convex Hull, Kirkpatrick and Seidel [9 described an algorithm 
to compute the Convex Hull of a set of n planar points in time within 0(n(l + log/i)) C 0(n\ogn), where 
h is the number of vertices in the Convex Hull. The algorithm relies on a variation of the divide-and- 
conquer paradigm, which they call the “Marriage-Before-Conquest” principle. For computing the upper hull, 
the algorithm finds a vertical line that divides the input point set into two approximately equal-size parts 
in liner time. Next, it determines the edge of the upper hull that intersects this line in linear time. It then 
eliminates the points that lie underneath this edge and finally applies the same procedure to the two sets of 
the remaining points on the left and right side of the vertical line. A similar algorithm computes the lower 
hull. Afshani et al. [l refined the complexity analysis of this algorithm to within 0(n( 1 + 7~L(ni, ..., n^))) C 
0(n(l+ log h)) C 0(n\ogn), where n\,... ,nh are the sizes of a partition of the input, such that every element 
of the partition is a singleton or can be enclosed by a triangle whose interior is completely below the upper 
hull of the set, and l~L(n\, ... ,nh) has the minimum possible value (minimum entropy of the distribution of 
the points into a certificate of the instance). This result takes advantage of the positions of the points i.e., 
the structure of the instance. 


2.2 Input-Order Based Results 


By “Input-Order Based Results” we mean algorithms taking advantage of the order of the input, for example, 
taking advantage of the order of the values in a sequence of numbers or of the order in which the points 
are given in a polygonal chain. We describe only a sampling of such results on the problem of Sorting 


permutations (see the survey from Moffat and Petersson 13 for more), revisit some results on the compu¬ 
tation of Convex Hulls in the plane and 3D space, and show that those results in the plane are actually 
only “Input-Order Based”. Those results for computing Convex Hulls in 3D space show that no algo¬ 
rithm can take advantage of the position of the points i.e., structure based, in order to compute Delaunay 
triangulations in the plane. 

Concerning the problem of Sorting, Knuth 


10 described an adaptive sorting algorithm that takes 


advantage of permutations formed by sorted blocks called runs, that is, subsequences of consecutive positions 
in the input with a positive gap between successive values, from beginning to end. He showed that the time 
complexity of this algorithm is within 0(n(l+ logp)) C 0(n\gn), where p is the number of runs in the 
permutation (e.g. (1,2,6,7,8,9, 3,4,5) is composed of 2 such sorted blocks (1, 2, 6, 7,8,9) and (3,4,5)). Barbay 
and Navarro [2 refined the algorithm described by Knuth jlO . They included in the analysis not only the 
number of runs but also their sizes, The main idea is to detect the runs first and then merge them pairwise, 
using a mergesort-like step. The detection of ascending runs can be done in linear time by a scanning process 
identifying the positions i in n such that n(i) > n(i + l). Merging the two shortest runs at each step further 
reduces the number of comparisons, making the running time of the merging process adaptive to the entropy 
of the sequence of the lengths of the runs. The merging process is then represented by a tree with the shape 
of a Huffman \7 tree, built from the distribution of the runs sizes. They extend this result to mix ascending 
and descending runs, showing that, if the permutation 7r is formed by p runs of sizes given by the vector 
(ri,... ,r p ), then n can be sorted in time within 0(n(l-\-H(r\,... ,r p ))) C 0(n(l+logp)) C O(nlogn). This 
result takes advantage of the order of the values in the input i.e., the input order. 

Considering the computation of the Convex Hull in the plane, Levcopoulos et al. E3 described a divide- 
and-conquer algorithm for computing the Convex Hull of a polygonal chain. The algorithm is based in the 
fact that the Convex Hull of a simple chain can be computed in linear time, and that deciding whether a 
given chain is simple can be done in linear time. They measured the complexity of this algorithm in terms 
of the minimum number of simple subchains n into which the chain can be cut. They showed that the 
time complexity of this algorithm is within 0(n(l+ log n)) C O(nlogn). We improve the analysis of this 
algorithm including not only the minimum number of simple subchain into which the polygonal chain can 


be partitioned but also their sizes (see Section 3.2). This result takes advantage of the order in which the 
points are given i.e., the input order. 

Given a set P of n planar points, a triangulation of P is a subdivision of the convex hull of P into 
triangles with vertex set the set P. Concerning the computation of the Convex Hull in 3D, a related 
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Fig. 1. A point set P. a) The convex hull of P, b) the Delaunay triangulation of P, and c) the Delaunay triangulation 
and the Voronoi diagram of P. 


concept is that of the Delaunay Triangulation DT(P) of a point set P in the plane, a triangulation 
where for every edge e there exists a disk C with the following properties: (i) the endpoints of edge e are on 
the boundary of C, and (ii) no other point of P is in the interior of C: it is named after Boris Delaunay for 
his work on this topic in 1934. An equivalent definition is such that no point in P is inside the circumcircle 
of any triangle of DT(P). Computing the Delaunay Triangulation is equivalent to computing its dual, 
called the VORONOI Diagram: each one can be constructed from the other in linear time jl5] (see Figure [7]). 
Computing the Delaunay Triangulation of a set of points in two dimensions reduces to computing the 
Convex Hull in three dimensions of the projections of those points on an hyperbolic plane. The projection 
of P onto the unit elliptic paraboloid z = x 2 + y 2 yields a point set P'. The Convex Hull CH(P') of 
P' contains every point of the set. The downward-facing facet of CH(P') are those whose normal vectors 
have a negative z- value. Projecting the edges of downward-facing facet in CH(P’) onto the plane yields the 
Delaunay triangulation of P. By showing tight bounds for input-order oblivious i.e., structure based, 
algorithms for computing Convex Hulls in three dimensions, Afshani et al. [l] indirectly proved that no 
planar Delaunay Triangulation algorithm can take advantage of the position of the points. 

Theorem 1 (Afshani et al. 0 )- Consider a set of n points in the plane. For any algorithm A computing 
the Delaunay triangulation in the algebraic decision tree model , A performs in time within /2(nlogn) on 
average on a random order of the points. This implies that there is an order of those points for which A 
performs in time within f2(nlogn). 

We describe in Sections [3] and [4] some algorithms taking advantage of the order of the input to compute 
Delaunay Triangulations and Voronoi Diagrams, among other desirable objects in computational 
geometry. 


3 Refined Analysis: Three Examples 

We show in this section how to refine the analysis of the principal step in the algorithm for multiplying 
polynomials using the Fast Fourier Transformation 16 , how to refine the analysis of the algorithm from 
Levcopoulos et al. 11 for the decomposition of a polygonal chain into simple sequences, and how to extend 


this analysis to the computation of VORONOI DIAGRAMS and Delaunay triangulations for another 
measure of difficulty based on monotone histograms. 


3.1 Polynomial Multiplication: Adaptivity to Zero-coefficients 

Given two polynomials A = (ao,..., a„_i) and B = (bo, ■ ■., b n - 1 ) described by their coefficients, the poly¬ 
nomial multiplication problem is to compute the coefficients of the polynomial C = A ■ B. The approach 







to multiplying polynomials using the Fast Fourier Transformation 16 can be divided into three steps: (if 
evaluate the polynomials A and B in 2n values (the (2n)th roots of the unity); (ii) evaluate C in these 2 n 
values by multiplying the evaluations of A and B ; and (iii) obtain the coefficients of C by interpolation 
using the values computed in the step (ii). The steps (i) and (iii) are accomplished by a divide-and-conquer 
algorithm for polynomial evaluation. We refine the analysis of the divide-and-conquer polynomial evaluation 
algorithm to take advantage of the number of zero-coefficients and of their relative positions in the vector of 
coefficients that describes the polynomial. 

Given a polynomial A, the polynomial evaluation algorithm defines two polynomials, A even and A a dd, 
that consist of the even-indexed and odd-indexed coefficients of A, respectively. Hence, A[x) = A even (x 2 ) + 
x x A a dd(x 2 ). If x is one of the (2n)th roots of the unity, then x 2 is one of the ?rth roots of the unity. In order 
to evaluate the polynomial A on each of the (2n)th roots of the unity, the recursive procedure divides A into 
A even and A a dd , evaluates A even and A 0 dd in the nth roots of the unity, and once these values are computed, 
evaluate A on each of the (2n)th roots of the unity using the formula A{x) = A even (x 2 ) + x x A a dd(x 2 ). 
The time complexity T(n) of this algorithm follows the recurrence T[n ) < T(^) + O(n), which yields a time 
complexity within 0(n log n). If at one step of the recursion call all the coefficients are zero, the algorithm 
finishes the computation at this branch (see Figure |2|. 


(4,0,0,3,2,0,0,1) 

(4,0,2,0) (0,3,0,1) 

(4,2) (0,0) (0,0) (3,1) 

(4) (2) (3) (1) 

Fig. 2. Vector formed by the coefficients of the polynomial A(x) = 4 + 3x 3 + 2x 4 + x 7 and the recursion tree of the 
divide-and-conquer evaluation algorithm. At each step, the algorithm divides the coefficients of A into A even and 
A 0 dd' 


Given a polynomial A = (do,..., a„_i), we define the equivalence relation E between the positions of the 
zero-coefficients in A: two positions of zeros p and q are equivalent if and only if (i) there exists k € Z + such 
that p = q (mod 2 fc ); (ii) all the positions r such that r = p (mod 2 k ) are zeros; and (iii) there exist a position 
t such that t = p (mod 2 k ~ 1 ) and the value in t is different from zero. The idea of this equivalence relation 
is to group the positions in classes such that the positions in the same class form a vector where all the 
coefficients are zero in a node of the recursion tree. Let £ and be the number of zeros in the vector formed 
by the coefficients of A and the number of equivalence classes defined by E, respectively. Let (m,... ,n v ) 
be the vector formed by the sizes of the equivalence classes defined by E. Then, l = C- The following 
theorem sets the refined analysis in function of ( and the vector (ni,..., n v ). 

Theorem 2. The complexity of the evaluation algorithm for polynomials of n coefficients, number of zero- 
coefficients f and vector (n \,..., n v ) formed by the sizes of the equivalence classes defined by E is within 
0{{n — Q log n + n i lo S ^ 0{n\ogn). 

In the following sections, we apply similar techniques to obtain optimal refined analysis for input-order 
adaptive algorithms computing Convex Hulls and Delaunay triangulations in the plane. 

3.2 Computing Convex Hulls: Adaptivity to Simple Subchains 

A polygonal chain is a curve specified by a sequence of points pi,p 2 , ■ ■ ■ ,p n . The curve itself consists of the 
line segments connecting the pairs of consecutive points. A polygonal chain C is simple if any two edges of 
C that are not adjacent are disjoint, or if the intersection point is a vertex of C; and any two adjacent edges 



ihare only their common vertex. Melkman 12 described an algorithm that computes the Convex Hull 
of a simple polygonal chain in linear time, and Chazelle [§] described an algorithm for testing whether a 
polygonal chain is simple in linear time. 

Levcopoulos et al. 11 combined these results to yield an adaptive divide-and-conquer algorithm for 
computing the Convex Hull of polygonal chains. The algorithm tests if the chain C is simple, using 
Chazelle '3i|’s algorithm: if the chain C is simple, the algorithm computes the Convex Hull of C in 
linear time, using Melkman 
subsequences C' and C' 


12 


s algorithm. Otherwise, if C is not simple, the algorithm cuts C into the 
whose sizes differ at most in one; recurses on each of them; and merges the 

They measured the complexity 


resulting CONVEX Hulls using Preparata and Shamos’s algorithm 15 


of this algorithm in terms of the minimum number of simple subchains n into which the chain C can be 
cut. Let t(n, k) be the worst-case time complexity taken by this algorithm for an input chain of n vertices 
that can be cut into n simple subchains. They showed that t(n, n) satisfies the following recursion relation: 
t(n,n) < t(|~f], Ki)+f(|_f J , K2), K1+K2 < k+1. The solution to this recursion gives t(n, k) € 0(ra(l+log/c)) C 
0(n\ogn). In the sequel, this algorithm will be named as Test-And-Divide. 

The Test-And-Divide algorithm partitions the input chain into simple subchains. If it was possible to 
partition the input chain into the minimum number k of simple subchains in linear time, then the same 
approach described by Barbay and Navarro [ 2 ] could be applied to obtain a refined analysis in function of 
0(n(l + ... ,n K ))) C 0(n(l+log«)) C O(nlogn). But, as far as we know, there does not exist any 

linear time algorithm to accomplish this task. 

In the recursion tree of the execution of the Test-And-Divide algorithm on input C of n points, every 
node represents a subchain of C. The cost of every node is linear in the size of the subchain that it represents. 
The simplicity test and the merge process are both linear in the number of points in the subchain. When 
this subchain is simple the node that represents this subchain is a leaf. Every time this algorithm discovers 
that the polygonal chain is simple, it executes a number of operations linear in the size of the chain and the 
corresponding node in the recursion tree becomes a leaf. 


Width Analysis: A Warm-up. Consider for illustration the particular case of a polygonal chain C of 
n = 2 m planar points such that C can be partitioned into the minimum number of simple subchains k = m 
of lengths 2 1 ,2 1 ,2 2 , 2 3 ..., 2 m_1 , respectively. Every time the algorithm cuts the current chain in half, the 
right subchain is simple and the recursive call is made only in the left subchain. Hence, the recursion tree 
of the algorithm on input C has only two nodes per level, one of which is a leaf (see Figure [3]). The overall 
running time of the algorithm on C is then within 0(n). 


0(f) 0(f) 

O(f) O(f) 

O(f) 


0 ( 1 ) 0 ( 1 ) 


Fig. 3. The recursion tree of A on C. Each node represents a recursive call. Noted in each node is the asymptotic 
complexities of the simplicity test and the merging process on the subchain that it represents. 


Definition 1 (Width). The width uj of the recursion tree in the execution of the Test-And-Divide algo 
rithm on input C is the maximum number of nodes at any level. 






Levcopoulos et al. 11 analyzed the complexity of this algorithm in the worst case over instances of fixed 


size n and n. The following lemma gives an alternate analysis in the worst case over instances of fixed size 
n and width ui. 


Lemma 1. Let to be the width of the recursion tree in the execution of the Test-And-Divide algorithm on 
input C of n planar points. The complexity of this algorithm on input C is within 0(um). 


Refined Analysis. Let (ii,... ,£ m ) be the vector formed by the sizes of the subchains represented by the 
m leaves of the recursion tree of the Test-And-Divide algorithm on input C (such that ^"= 1 ^ = n )• The 
number of operations “saved” by the algorithm every time it discovers a leaf of size ti is within 12 {ii log if) 
(the cost of the subtree of the perfect binary tree rooted in a node of size if) minus 0(£f) (the operations in 
the leaf are not saved). The time complexity T(C) of this algorithm on input C is within O(nlogn) (the cost 
of the perfect binary tree) minus L2(J2"f 1 £ i logii — if) (the number of operations saved by the algorithm). 
So, T{C) C 0(n log n — X)™ 1 (£ t log it-if)) = 0(n(l + TL(t\,... ,£ m ))) C 0(um) D O(nlogm) C O(nlogn). 

The following lemma summarizes this finer analysis of the Test-And-Divide algorithm in function of the 
width and the number of leaves in the recursion tree. 

Lemma 2. Let £ m ) be the vector formed by the sizes of the subchains represented by the m leaves of 

the recursion tree in the execution of the algorithm A on input C of n = £i planar points. Let w be the 
maximum width of the recursion tree. The time complexity of A on C is within 0{n{ 1 + TL{i i,... ,£ m ))) C 
0{um) (~l 0{n log in) C 0(n log n). 

Is there a relationship between the vector {i\,... ■ £ m ) formed by the sizes of the subchains represented 
by the leaves of the recursion tree in the execution of the Test-And-Divide algorithm from on input C and 
the vector (m,..., n K ) formed by the sizes of a partition of C into k simple subchains? 

For a given polygonal chain, there can be several partitions into simple subchains of minimum size k 
for it. The Levcopoulos et al.’s analysis is in the worst case over instances for n and n fixed. We describe 
below a refined analysis which takes into account the relative imbalance between the sizes of the subchains. 
The idea behind the refinement is to bound the number of operations that the algorithm executes for every 
simple subchain. This analysis makes it possible to identify families of instances where the complexity of the 
algorithm is linear even though the number of simple subchains into which the chain is split is logarithmic. 

Theorem 3. Let («i,..., n K ) be the vector formed by the sizes of the subchains of any partition II of the 
chain C into the minimum number n of simple subchains. The time complexity of the Test-And-Divide 
algorithm on input C is within 0{n{ 1 + TL(n \,..., n K ))) C 0(n( 1+ log k)) C 0(n logro), which is worst-case 
optimal in the comparison model over instances of n points that can be partitioned into n simple subchains 
of sizes (ni,..., n K ). 

Proof. Fix the subchain Cj of size Ui in 77. In the worst case, the algorithm considers the n,; points of Ci for 
the simplicity test, and the merging process in all the levels of the recursion tree from the first level to the 
level [log —] +1, because the sizes of the subchains in these levels are greater than n,. In the next level, one 
of the nodes t of the recursion tree fits completely inside c* and therefore it becomes a leaf. Hence, at least 
^ points from Ci are dismissed for the following iterations. The remaining points of Ci are in the left or the 
right ends of subchains represented by nodes in the same level of t in the recursion tree. In all of the following 
levels, the number of operations of the algorithm involving points from Cj can be bounded by the size of the 
subchains in those levels. So, the sum of the number of these operations in these levels is within 0{nf). As 
a result, the number of operations of the algorithm involving points from Ci is within 0(ni log — + nf). In 
total, the time complexity of the algorithm is within 0(n + l n i^°S )yr) = 0(n( 1 + ... ,n K ))) C 

0(n(l+logre)) C O(nlogn). 

We prove the optimality of this complexity by giving a tight lower bound. Barbay and Navarro [2| showed 
a lower bound of I2(n(l + 77(r 1; ... ,r p ))) in the comparison model for SORTING a sequence of n numbers, in 
the worst case over instances covered by p runs (increasing or decreasing) of lengths n,..., r p , respectively, 
summing to n. The SORTING problem can be reduced in linear time to the problem of computing the 



fcoNVEX Hulls of a chain of n planar points that can be cut into p simple subchains of lengths ri,...,r p , 
respectively. For each real number r, this is done by producing a point with (x, y)-coordinates (r, r 2 ). The 
p runs (alternating increasing and decreasing) are transformed into p simple subchains of the same lengths. 
The sorted sequence of the numbers can be obtained from the Convex Hull of the points in linear time. □ 


3.3 Computing Delaunay Triangulations: Adaptivity to Monotone Histograms 

We propose a new input-order based algorithm for computing the Delaunay TRIANGULATION and the 
VORONOI DIAGRAM of a set of n planar points. A monotone histogram is a sequence of points sorted with 
respect to two orthogonal directions. The algorithm takes advantage of the minimal number p of monotone 
histograms and their sizes into which a polygonal chain can be cut. A monotone histogram is also a simple 
polygonal chain. Therefore, an algorithm for computing Convex Hulls adaptive to the decomposition of 
the input into monotone histograms is obtained as a corollary of Theorem [3j We extend the refined analysis 
to the computation of Delaunay triangulations and Voronoi diagrams adaptive to the decomposition 
of the input into monotone histograms. 

Djidjev and Lingas [4 described an algorithm which, given a monotone histogram, computes the VORONOI 
DIAGRAM (and hence the Delaunay triangulation) of the input sequence in linear time. This algorithm 
suggests a way to partition the input into subsequences such that the Delaunay triangulation of each 
subsequence can be computed in linear time in its length. The partitioning algorithm cuts the sequence into 
p monotone histograms with respect to two orthogonal directions d\ and c? 2 - The first two points of the 
subsequence determine the ordering defined by the combination of ascending and descending with respect 
to d\ and e? 2 • Given a sequence of points and two orthogonal directions, it is possible to test whether the 
sequence is a monotone histograms with respect to these two directions in linear time. 


Binary Merge of Voronoi Diagrams. Kirkpatrick [8] described a linear time algorithm for the merging 
of two arbitrary Voronoi DIAGRAMS. Given the VORONOI DIAGRAMS of two disjoint point sets P and Q , 
the algorithm finds the VORONOI DIAGRAM of P U Q in time within 0(|P| + |Q|). The plane is partitioned 
into points closer to P, points closer to Q , and points equidistant from P and Q. The points equidistant 
from P and Q are defined as the contour separating P and Q. The contour is composed of straight line 
segments: it is formed from the edges of the VORONOI DIAGRAM of P U Q that separates the points in P 
from the points in Q. Inside the region of points closer to P (resp. Q) the VORONOI DIAGRAM of PU Q and 
the VORONOI diagram of P (resp. Q) are identical. Thus, the merging of two VORONOI DIAGRAMS can be 
seen as the process of cutting the VORONOI DIAGRAMS of P and Q along the contour. 

This leads to a divide-and-conquer algorithm for constructing the VORONOI DIAGRAM of a set of n points 
and hence for computing the Delaunay TRIANGULATION in time within O(nlogn). This time complexity 
of 0(n log n) is asymptotically optimal in the comparison model in the worst case over instances composed 
of n points. Shamos and Hoey 17 showed that the construction of any triangulation over n points requires 
Q(n log n). as Sorting can be reduced to computing the triangulation of n + 1 points, which yields an 
asymptotic computational lower bound of I2(nlogn) in the worst case over sets of n planar points, in the 
comparison model. 


Multiary Merge. Given p Delaunay triangulations of sizes iq, ... , 1 V, respectively, to be merged, 
we make a sequence of binary merges, reducing at each step the number of Delaunay triangulations 
to be merged by 1. The merging process can be represented by a binary tree where the internal nodes are 
the merged Delaunay triangulations, and the leaves are the original p Delaunay triangulations. 
Merging the two Delaunay triangulations of minimum sizes at each step, similar as the Huffman code 
algorithm |7 , further improves the merging process, which takes advantage of the potential disequilibrium 
in the distribution of the points between the p Delaunay triangulations. We can apply the Huffman [T] 
algorithm to the vector (tq,..., v p ), thus obtaining a Huffman-shaped tree representing the merging process. 



Lemma 3 (Multiary Merge). Given p Delaunay triangulations of respective sizes {v\,...,v p ) sun$ 
ming to n = 1 v i> there is an algorithm computing the Delaunay triangulation of the n points in 

time within 0(n( 1 + LL(v 1 ,..., v p ))) C 0(n(l + log/i)) C O(nlogn). 

Proof. The algorithm follows the same steps as the algorithm suggested by Huffman 7 on a set of p messages 
of probabilities {D/n}ie[i..p] : 

1. Initialize a heap H with the p Voronoi DIAGRAMS, indexed by their size; 

2. While H contains more than one VORONOI DIAGRAM 

— extract the two smallest VORONOI DIAGRAMS from H, of respective sizes n\ and ri 2 , 

— merge them into a VORONOI DIAGRAM T of size ni + ri 2 , and 

— insert T in H. 

This algorithm executes in time within 0(n(l + TL(v \,..., v p ))). The merging process is then represented 
by a tree with the shape of a Huffman |7j tree. Consider the i-th Voronoi DIAGRAM of the input Vi £ [1..//]: 
let Ci be the binary string describing the path leading from the root to the corresponding leaf, and f the 
length of this path. The sum of the computational costs of the binary merges is the sum of the sizes of the 
Voronoi diagram computed. The i-th Voronoi diagram contributes a cost within 0(vi) to k levels, 
which has a sum within 0(X]f=i hvi). Consider the binary tree where the p initial VORONOI DIAGRAMS are 
leaves and the p — 1 computed VORONOI DIAGRAMS are internal nodes: 

— The set of binary strings {c \,..., c p } is a prefix free code, i.e. no code is prefix of another root-to-leaf 
path, simply because they are paths in a tree. 

— The lengths of those codes minimize )T^ =1 hXi as a property of Huffman [ 7 ] codes. 

By the optimality of Huffman codes, this complexity is within a linear term of the entropy of the distribution 
(ui,..., v p ), i.e. X]f=i li r i £ 0(n(l + TL(v\,..., v p ))). This yields the final time complexity, within 0(n( 1 + 
'H{V Vfj))). □ 

The combination of the partitioning algorithm and the merging process yields an optimal algorithm 
computing these structures adaptive to the decomposition of the input into monotone histograms. 

Theorem 4. Let d\ and d 2 be two perpendicular directions. Let S be a sequence ofn planar points. Let p and 
{v \,..., vP) be the minimum number of monotone histograms with respect to g?i and d 2 and the sizes of these 
monotone histograms, respectively, in which S can be cut. The Delaunay triangulation and the Voronoi 
diagram of S can be computed in time within 0(n(l + TL(vi,... ,v p ))) C 0(n{l+ logp,)) C O(n logn), which 
is worst-case optimal in the comparison model over instances ofn points that can be partitioned into monotone 
histograms of sizes (v \,..., v^). 

Proof. The combination of the partitioning algorithm and the merging process yields an algorithm computing 
these structures within this time. In order to provide a lower bound, we use again the result of D(n(l + 
TL(ri,... ,r p ))) for Sorting a sequence of n numbers, in the worst case over instances covered by p runs 
of lengths r\,... ,r p summing n in the comparison model, demonstrated by Barbay and Navarro [2j. This 
problem can be reduced in linear time to the problem of computing the Delaunay triangulation of 
a sequence of n planar points covered by p monotone histograms of lengths r\,...,r p with respect to the 
coordinates axes. The p runs are transformed into p monotone histograms of the same lengths, using points 
on the parabola, in linear time. The sorted sequence of the numbers can be obtained from the Delaunay 
triangulation of the points in linear time. □ 

Such examples of complete refinements are not the rule: in the following section we describe some examples 
where such refinements are problematic or impossible. 

4 Partial Refinements 

We describe a new partitioning algorithm for a sequence of points in Sections |4.1| to |4.3| such that the 
Delaunay triangulation of each subsequence can be computed in linear time in its length, and discuss 
alternate partitions in Section |4.4| Each different partitioning algorithm, in combination with the merging 
process described in Section |3.3[ yields a different algorithm adaptive to the input order. 


49l Incremental Construction 


Many Computational Geometry algorithms use tests known as the orientation and incircle tests 18 . The 


orientation test determines whether a point lies to the left of, to the right of, or on a line or plane defined 
by other points. The incircle test determines whether a point lies inside, outside, or on a circle defined by 
other points. 

Green and Sibson [5] proposed the first incremental algorithm for computing the Delaunay trian- 
GULATION of a point set which finds the triangle containing each new point, and updates the diagram by 
correcting the edges violating the circumcircle condition —an operation named flipping. The algorithm adds 
points to the structure one by one. For each point, it performs two basic steps: 


1 . 


2 . 


The algorithm finds the triangle containing the new point using the structure as a guide to the relative 
position of the points. A greedy approach for locating the point is to start at an edge in the structure 
and to walk across adjacent edges in the direction of the new point until the correct triangle is found. 
Orientation tests 18 are performed on each edge of such a path to see whether the new point lies on 


the correct side of that edge. We call the operations involved in this walk navigation operations. 

It then updates the structure adding three new edges from the point inserted to the vertices of the 
triangle containing the point and flips all invalid edges resulting from the insertion. Note that flipping 
an edge can make another edge invalid, but that each edge is flipped at most once, so each insertion can 
trigger at most a linear number of flips. 


The time complexity of this incremental algorithm for computing the Delaunay TRIANGULATION is within 
0(dn ) C 0(n 2 ), where d £ [l..n] is an upper bound on the amount of operations required to insert each point 
and to correct the Delaunay triangulation for the instance being considered. We show the problems 
that arise when we try to adapt this algorithm in order to obtain a new one whose time complexity is within 
O(nlogd) C 0{n log n). 


4.2 Adaptivity to D-linear Runs 

Consider an incremental algorithm G computing the Delaunay triangulation of n planar points in time 
within 0(n) in the best case; and an algorithm M merging two Delaunay triangulations of sizes summing 
to n in time within 0(n) in the worst case. Given a sequence S of n distinct planar points and a constant k, the 
following naive algorithm computes its Delaunay triangulation in time within 0{n ) in the best case, in 
time within 0(n log n) in the worst case, and in time within 0(n(l+'H(ri, ..., r p ))) C O(nlogp) C O(nlogn) 
in the general case, where p £ [1 ..n\ and ri,...,r p (such that n = r *) measure the difficulty of 

partitioning the instance into “easy” subinstances: 

1. Run the incremental algorithm G on S until, for the i-th point p , it performs either more than k operations 
to locate the point on the triangulation, or more than k flip operations on it. 

2. Store the Delaunay triangulation of the i — 1 first points, and restart the greedy incremental 
algorithm G on the sequence formed by p and its n — i successors in the input sequence S, until no points 
are left in S. 

3. Let p £ [l..n] and n,...,r p (such that n = J2i=i r i) be the number of Delaunay triangulations 
computed in this way and the sizes of these Delaunay triangulations, respectively. Merge the p De¬ 
launay triangulations, in overall time within 0(n(l+'H(ri,... ,r p ))) C 0(n(l+logp)) C O(nlogn). 


Greedy Partitioning. The algorithm described above suggests a way to partition the input into subse¬ 
quences such that the Delaunay triangulation of each subsequence can be computed in linear time in 
its length. 

We use the algorithm described to define “easy” sequences of points for the computation of the Delaunay 
TRIANGULATION: 
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Fig. 4. The sequences A, B and A'. 


Definition 2 (D-linear). Consider a sequence S ofn distinct planar points. Given an integer value k > 0, S 
is /c-D-linear (for “Delaunay linear”) if G performs at most k location and flip operations for each point while 
computing the Delaunay triangulation of S. If S is k-D-linear and k G 0(1) is a constant independent 
ofn, we say that S is D-linear. 

Such a simple definition yields a simple partitioning of the input sequence into subsequences of consecutive 
positions such that the Delaunay triangulation of each subsequence can be computed in linear time. 

Definition 3 (D-linear Run). A D-linear run in a sequence S of points is a D-linear subsequence formed 
by consecutive points in S. 

4.3 Non Optimality of D-linear Runs Partitioning 

We describe a family of sequences S, such that for all positive integer n, there is a sequence S n of n 
points in S where the greedy partitioning algorithm yields three D-linear runs when S n can be optimally 
partitioned into just two D-linear runs. 

The Delaunay triangulation of m > 2 points on the parabola y = x 2 is such that the leftmost point 
is adjacent to all the other m — 1 points, and the left-to-right order is a D-linear run. The proof of the 
following lemma is based on this fact. 

Lemma 4. For all positive integers n, there exists a sequence ofn planar points where the greedy partitioning 
algorithm yields three D-linear runs, whereas the optimal partition of this sequence has only two D-linear 
runs. 

Proof. Let k be the threshold used in the definition of D-linear sequences. Let A be a set of m > k points in 
the parabola y = x 2 , denoted pi, ■ ■ ■ ,p m from left to right. Let C be the circumcircle of the triangle Pi,P 2 ,Ps- 
Let u , v, w be three points such that the triangle Auvw with vertex set { u , v, w} is small enough with respect 
to the convex hull of A , and located inside C. Let A! = {p[,... , p' rn } be a scaled copy of A, located inside 
Auvw. Let B be the k + 3 vertex set of k + 1 decreasing area and disjoint triangles such as every two 
consecutive triangles share an edge (see Fig. [4]). Suppose that B is ordered such that the sequence of points 
is a D-linear run. The points u, v and w are the last three points of B. The triangle formed by the first three 
points of B contains the point p\ of A. Let P = AcA'UB be ordered as (B,pi,..., p m , p\ , ■ ■ ■ ,p' m )- Applying 
the greedy partitioning algorithm to P gives the three D-linear runs (B), (pi,... ,p m ), and {p\,... ,p' m ) since: 
(1) adding p\ to the Delaunay triangulation of B requires traversing k triangles to locate p±; and (2) 
adding p\ to the Delaunay triangulation of A requires creating m triangles (i.e. p\ is adjacent to every 
point in A). However, (B \ {u,v,w}) is a D-linear run and (u,v,w,pi,...,p m ■ p \, ■ ■ ■ ,p' m ) is another one. 
Note that the triangle Auvw blocks the points in A! of being adjacent to any point in A. □ 

This family of sequences S can be generalized to show a gap of 0{n) between the number of D-linear 
runs yielded by the greedy partitioning algorithm in S n and the minimum number of D-linear runs in which 
S n can be partitioned. 






&mma 5. For all positive integer n, there exists a sequence of n planar points where the greedy partitioning 
algorithm yields p £ 0(n ) D-linear runs, and the optimal partition of this sequence has only 2 D-linear runs. 

Proof. Let A be a set of m > k consecutive points in the parabola y = x 2 denoted Pi, ■ ■ ■ ,p m from left to 
right. Let C be the circumcircle of the triangle p\,p 2 ,P 3 - Let u,v,w be three points such that the triangle 
Auvw with vertex set { u , v, w} is small enough with respect to the convex hull of A , and located inside C. 
Let A' be a D-linear run sequence of l points, of increasing y-coordinate and decreasing x-coordinate , inside C 
denoted p\,...,p\. During the incremental construction of the Delaunay triangulation of the sequence 
AA! the first point of A! forces a flipping of m edges in A (i.e. m new edges are created connecting p[ with 
all points in A). Every next point of A' forces a flipping of the m edges connecting points in A with points 
in A'. Let q[ and q 2 be two consecutive points from the sequence A' and let q^,..., q k +i be k + 1 consecutive 
points from A. The sequence Q = (q[, q ±,..., qk+i,q' 2 ) is not a D-linear run because q 2 forces a flipping of the 
k + 1 edges connecting q[ with all the points in q\,q 2 ,..., q k +i- Let B be the [k + 3) vertex sequence of k + 1 
decreasing area and disjoint triangles so that every two consecutive triangles share an edge. The points u, v 
and w are the last three points of B. The triangle formed by the first three points of B contains inside the 
point pi of A. B is ordered such that the sequence of vertices is a D-linear run. Let P = dUA'UB be ordered 
as ( B,pi ,... ,Pk+i,Pi,Pk+ 2 , ■ ■ ■ ,P 2 k+ 3 ,P 2 > ■ ■ ■ ? Pm): that is, first the points of B , then alternate k + 1 points 
of A with 1 point of A', resulting subsequences similar to Q. The greedy partition algorithm divides P into 
min(7, D-linear runs: (B), (pi,... ,pk+i) and every subsequence starting with 1 point from A’ followed 

by k + 1 points from A. Since: (i) adding p\ to the Delaunay triangulation of B requires traversing k 
triangles to locate pi; (ii) adding p[ to the Delaunay triangulation of p±.. -Pk+i requires creating k 
triangles (i.e. p[ is adjacent to every point in p x ...p k+1 ); and (iii) adding p' i+1 to ( p[,Pik+j, ■ ■ ■ ,P(i+i)k+j+i) 
produces a sequence similar to Q , and hence is not a D-linear run. However, (B \ {zt, v, w}) is a D-linear run 
and (zt, v, w,pi,... ,Pk,Pi,Pk+i>Pk+ 2 , ■ ■ ■ ,P 2 k,P 2 i ■ ■ ■ iPm) is another one. Note that the triangle Auvw blocks 
the points in A' of being adjacent to any vertex in A. It is possible to build a sequence l = -jffj = n f^_ 2 3 
(1 point in A' for each sequence of k + 1 points in A) such that the number of D-linear runs yielded by the 
greedy partitioning algorithm is n ~f^_ 2 3 + 1 and the optimal partition has only 2 D-linear runs. □ 

4.4 Other Partitioning Schemes 

Since the greedy partitioning algorithm does not yield an optimal partition, we explore more sophisticated 
partition techniques and show that they suffer similar problems. 


We adapt the Levcopoulos et al. 11 partition technique to cut a sequence of planar points into D-linear 


runs. The test and divide partitioning is another partition technique based on the incremental algorithm G for 
the computation of the Delaunay triangulation (seen in Sect. 4.11. Given a sequence S = (pi,. ■ ■ ,p n ) 
of n planar points, the test and divide partitioning first tests whether the sequence S' is a D-linear run. 
In such a case, it identifies the sequence S as a D-linear run. If not, it partitions the sequence S into 
S' = (pi,... ,P[ n / 2 j) and S" = {p\ n / 2 \+i, ■ ■ -iPn) and recursively the same procedure is applied to S' and 
S". While the greedy partitioning has a linear time complexity, this partitioning has a time complexity 
within 0{n log p) C 0(n log n) where p is the number of D-linear runs identified by the partitioning process. 

The following lemma shows that the test and divide partitioning suffers the same problems as the greedy 
partitioning. 


Lemma 6. For all positive integer n, there exists a sequence of n planar points where the test and divide 
partitioning yields p £ 0(n) D-linear runs and the optimal partition of this sequence has only 2 D-linear 
runs. 


Proof. Consider the sequence S used in the proof of Lemma[5]where the number of points in B is changed to 
|_n/2j instead of k + 3 then the test and divide partitioning will cut S in the last point of B. B is a D-linear 
run, but the bisection partitioning will cut the rest of S in 2 ^ k + 2 ) D-linear runs. □ 
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The greedy partitioning algorithm adds the point pi to the current run whether the number of location 
and update operations of the incremental algorithm G in pi is at most a threshold k. But it is possible that 
in some points the number r of operations would be much lower than k. The amortized greedy partitioning 
algorithm makes use of the k — r remaining operations in the subsequent points, similar to the accounting 
method for amortized analysis. The number k — r of remaining operations are credited to be used later, so 
that the point pi will be added to the current run if the number of navigation and update operations of G 
in pi is less than k plus this credit. 

We use this algorithm to give another definition of “easy” sequences of points for the computation of the 
Delaunay triangulation: 

Definition 4 (Amortized D-linear). Consider a sequence S = (p± : ... ,p n ) of n distinct planar points. 
Given an integer value k > 0, S is Amortized fc-D-linear if for each point pi £ S, the incremental algorithm 
G performs at most k * i navigation and flip operations while computing the Delaunay triangulation of 
the sequence Pi,... ,Pi- If S is Amortized k-D-linear and k £ 0(1) is a constant independent of n, we say 
that S is Amortized D-linear. 

Such definition yields a partitioning of the input sequence into subsequence of consecutive positions: 

Definition 5 (Amortized D-linear Run). An Amortized D-linear run in a sequence S of points is an 
Amortized D-linear subsequence formed by consecutive points in S. 

We define the optimal partition into Amortized D-Linear runs as the partition with the minimum number 
of Amortized D-linear runs. 


Lemma 7. For all positive integer n, there exists a sequence of n planar points where the amortized greedy 
partitioning yields p £ 0(n) Amortized D-linear runs and the optimal partition has only 2 Amortized D-linear 
runs. 


Proof. The construction of the sequence follows the same scheme of previous constructions (see the proof 
of Lemma [5]). The sequence B is such that the incremental algorithm G performs k navigation and flip 
operations in almost every point in B, so that B is still a run but the credit is almost zero at the end of B. 
Again, the last 3 points of B form a triangle that blocks the points in A' of being adjacent to any vertex in 
A. The sequence S alternates points in A with points in A! such that each new point in A' is adjacent to 
every point in A, thus when the numbers of points in A is large enough to exceed the credit, the Amortized 
D-linear run is cut. This partition algorithm yields close to n/k 2 Amortized D-linear runs. This sequence can 
be partitioned optimally in just 2 Amortized D-linear runs. Note that the triangle Auvw blocks the points 
in A' of being adjacent to any point in A. □ 


Amortized Test and Divide Partitioning. The amortized test and divide partitioning is a combination 
of the test and divide partitioning and the amortized greedy partitioning. Given a sequence S = (pi,... ,p n ) 
of n planar points we first test whether the sequence S is Amortized D-linear. In such a case, we identify the 
sequence S as an Amortized D-linear run. If not, we partition the sequence S into S' = (pi,... ,py n / 2 j) an d 
S" = {p\n/ 2 j+i> • • • iPnj and recursively the same procedure is applied to S' and S". 

Lemma 8. For all positive integer n, there exists a sequence of n planar points where the amortized test and 
divide partitioning yields p £ 0(n) Amortized D-linear runs and the optimal partition has only 2 Amortized 
D-linear runs. 


^roof. The construction of the sequence follows the same scheme of previous constructions (see the proof 
of Lemma [5]). The sequence B is such that the incremental algorithm G performs k navigation and flip 
operations in almost every point in B, so that B is still a run but the credit is almost zero at the end of B. 
Again, the last 3 points of B form a triangle that blocks the points in A' of being adjacent to any vertex in 
A. The sequence S alternates points in A with points in A! such that each new point in A! is adjacent to 
every point in A, thus when the numbers of points in A is large enough to exceed the credit, the Amortized 
D-linear run is cut. This partition algorithm yields close to n/k 2 Amortized D-linear runs. This sequence can 
be partitioned optimally in just 2 Amortized D-linear runs. Note that the triangle Auvw blocks the points 
in A! of being adjacent to any point in A. □ 


5 Discussion 


We describe one technique to refine the analysis of the divide-and-conquer polynomial evaluation algorithm, 
and two techniques to refine the analysis of divide-and-conquer input-order adaptive algorithms for computing 
the Convex Hull, Voronoi Diagram and Delaunay Triangulation of n planar points in time within 
0(n(l + 7-L(ni,... ,rik))) C 0(n(l+log/c)) C O(nlogn) for some parameters n\,...,nk summing to n, 
where H(ni ,..., rik) = Xu=i 7T 1°S rt measures the “difficulty” of the instance. We show the optimality of 
the algorithms for computing Convex Hulls, Voronoi Diagrams and Delaunay Triangulations by 
providing lower bounds that match the refined analyses. The later results improve on those from Levcopoulos 
et al. 11 in that this analysis takes advantage of the sizes of the simple polygonal chains into which the 


input polygonal chain can be partitioned. 

There is still work to be done in the directions initiated by our work. One of these techniques is based 
on the partitioning of a sequence of n planar points into subsequences, whose corresponding VORONOI DI¬ 
AGRAMS and Delaunay triangulations can be computed quickly, and a technique to efficiently merge 
those structures into a single one, inspired by Huffman |7| codes. Each combination of partition and merging 
techniques yields adaptive algorithms to the input order for computing VORONOI Diagrams and Delau¬ 
nay Triangulations. Those results show that the Delaunay Triangulation can be computed in time 
within o(nlogn) on some classes of instances, yet we described various partitioning algorithms where the 
time complexity achieved is sub-optimal. The next step is to find better partitioning algorithms taking full 
advantage of the input order. 

We classify the various refined analysis between those focusing on the structure of the instance versus 
those focusing on the order in which the input is given. We have been analyzing an algorithm for the Union 
Set problem that can be adapted for solving the Sorting problem taking advantage of the structure of the 
instance and the order of the input in synergy. We will continue studying potential synergistic solutions to the 
Sorting problem and their extension to the Convex Hull and Delaunay triangulation computation. 
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